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We study the transfer of cosmological perturbations through a nonsingular cosmological bounce 
in a special model in which the parameters of the bounce and the equation of state of matter are 
chosen such as to allow for an exact calculation of the evolution of the fluctuations. We find that 
the growing mode of the metric fluctuations in the contracting phase goes over into the growing 
mode in the expanding phase, a result which is different from what is obtained in analyses in which 
fluctuations are matched at a singular hypersurface. Consequences for Ekpyrotic cosmology are 
discussed in a limit when the equation of state of a fluid becomes large. 
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I. INTRODUCTION 

It has long been appreciated that Standard Big Bang 
cosmology cannot be the complete theory of the early 
universe since it is plagued by an initial singularity [ij. 
The singularity theorems have also been extended to ap- 
ply to scalar field-driven inflationary cosmology 01 . It is 
possible that the initial singularity is resolved by quan- 
tum gravity. Another possibility, however, is that cor- 
rection terms to the gravitational Lagrangian which de- 
scribes the dynamics of space-time leads to a bouncing 
universe. In a bouncing cosmology, the Hubble radius 
H~^(t) (where H is the Hubble expansion rate and t is 
physical time) decreases faster than the physical wave- 
length of fluctuations (which have a constant wavelength 
in comoving coordinates). Hence, it is conceivable that 
processes acting in the contracting phase can lead to a 
non-inflationary mechanism for the origin of structure in 
the universe. Recent proposals which involve bouncing 
cosmologies are the Pre-Big-Bang (3, Ekpyrotic and 
higher derivative modification of an Einstein gravity Q 
scenarios. 

The transfer of metric fluctuations through the bounce 
from the contracting to the expanding phase has been an 
outstanding problem in the discussions of Pre-Big-Bang 
and Ekpyrotic cosmology. Often, the evolution of the 
background is modeled by a contracting phase modeled 
by a solution to the Einstein equations matched to an ex- 
panding phase of Standard Big Bang cosmology through 
an instantaneous and often singular transition along a 
space-like hypersurface. One proposal has been to use 
the analog of the Israel matching conditions (matching 
conditions [6] which describe the merger of two solutions 
of the Einstein equations along a time-like hypersurface) . 
These equations were discussed in 0, Q . 

If we consider Einstein gravity, then in Fourier space 
the space of solutions for fixed comoving wavenumber 
k is two-dimensional (see e.g. @ for a comprehensive 
review of the theory of cosmological perturbations and 
for a pedagogical overview). We eliminate the coor- 
dinate ambiguities in the description of the fluctuations 



by working in a specific coordinate system, namely lon- 
gitudinal gauge, in which the metric is given by (in the 
absence of anisotropic stress) 

ds^ = a{Tf [(1 + 2$)dr2 - (1 - 2$)dx2] , (1) 

where t is conformal time related to physical time t via 
dt = adr, and the relativistic potential $(x, t) is the field 
describing the fluctuations. 

On super- Hubble scales and in an expanding universe, 
the dominant mode of $ is constant in time if the equa- 
tion of state of the cosmological background is constant. 
We will call this mode the D-mode. The second funda- 
mental solution, the S-mode, is decreasing in time. In a 
contracting phase, the D-mode is sub-dominant, and the 
S-mode is growing and hence dominant. 

Both in the context of Pre-Big-Bang scenario [ll| and 
in the original effective field theory description of the 
Ekpyrotic cosmology it was found that if the fiuctu- 
ations are matched across a hypersurface of fixed matter 
field value (the hypersurfaces singled out in the scenar- 
ios of d, 01), then the growing mode in the contracting 
universe couples almost exclusively to the subdominant 
mode in the expanding phase (the coupling to the dom- 
inant mode is suppressed by /c^). In the Ekpyrotic sce- 
nario [3] a scale-invariant spectrum of the dominant mode 
of cosmological fiuctuations is generated in the contract- 
ing phase [l3|, whereas the decaying mode has a n = 3 
spectrum ^ . Due to the suppression of the coupling, only 
a n = 3 spectrum in the expanding phase is induced ^ 



We are employing the standard notation in which the index of 
the spectrum of scalar metric fiuctuations is denoted by n — 1, 
with n = 1 standing for a scale-invariant spectrum. 
In the case of the Ekpyrotic scenario, the underlying physics is 
higher-dimensional, and the reduction of the analysis to a four 
space-time dimensional effective field theory without adding en- 
tropy modes does not correctly model the full physics of the 
fluctuations. An analysis of the transfer of metric fluctuations in 
five space-time dimensions IT3l (see also llSlI for a similar anal- 
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As initially pointed out in [T^, [13 , the result for the 
spectrum of fluctuations in the expanding phase depends 
quite sensitively on the details of the matching, and it 
is in fact not clear whether the matching prescription of 
0, Q can be applied at all due to the inconsistency of 
the matching of the background. It is thus of great inter- 
est to study how the cosmological fluctuations propagate 
through a nonsingular cosmological bounce. A first anal- 
ysis of the evolution of fluctuations through a specific 
nonsingular bounce obtained by using a higher deriva- 
tive gravity action was made in (see also [l^), show- 
ing that initial scale-invariant fluctuations do not pass 
through the bounce, thus confirming the results of [l2l |. 
On the other hand, bounces obtained by adding spatial 
curvature and matter with wrong-sign kinetic terms, but 
with the standard gravitational action, have been studied 
in [TgI. [pol. [2H ■ with differing results^. Whereas the anal- 
ysis of [2l| yielded no transfer of the growing mode in the 
contracting phase to the dominant mode in the expand- 
ing period, the analysis of ^16] showed that the transfer 
matrix which links the two fluctuation modes in the con- 
tracting phase to those in the expanding phase depends 
quite sensitively on the details of the bounce, and that it 
is possible that the coupling of the dominant mode in the 
contracting phase to that in the expanding phase is not 
suppressed. A similar conclusion was reached recently in 
a study of a single field model without spatial curvature 
in the context of a theory with modified kinetic term (23j 
^. In the context of a bouncing mirage cosmology it has 
also been shown very recently [23| that the spectral index 
of the dominant mode of $ does not change. 

The approach that we take instead, is to assume that 
the new physics which is responsible for resolving the 
singularity does not effect the evolution of the pertur- 
bations. The advantage of such an approach is that all 
the information of the new physics is now encoded in the 
background behaviour of the Hubble rate, and in par- 
ticular we can track the evolution of the perturbations 
unambiguously through the bounce, without having to 
implement any matching conditions. For some special 
cases one even has exact analytic expressions which fa- 
cilitates our understanding. Clearly, such an assumption 
cannot be valid for arbitrary new physics, but our main 
motivation for considering it comes from recent progress 
that has been made in resolving singularities involving 
"non-local" physics. In these cases one can argue that 
the new physics does not effect the "local" perturbation 
equation but only impacts the "global" or background 



ysis in a non-singular cosmological background) shows that in 
the higher-dimensional setup, an initial scale-invariant spectrum 
of cosmological perturbations in the contracting phase goes over 
into a scale-invariant spectrum in the expanding phase. 

^ Note that there is a technical problem with some of the analyses 
in these references, as discussed in 22]. 

4 See also ^ for an analysis of how fluctuations arise in a bounc- 
ing cosmology via quantum cosmological methods. 



evolution of the universe. For instance, in [26j a cosmo- 
logical BCS theory was considered where fermions can 
form Cooper pairs and the negative gap (binding) energy 
can mediate a bounce. Now the gap energy depends on 
the chemical potential of the system which in turn de- 
pends on the density of states. As is well known, the 
density of states in a given system depends only on the 
total volume of the system. Thus although the Gap en- 
ergy contributes to the Hubble equation, it is completely 
ignorant of local fluctuations of the volume (metric) and 
therefore does not alter the perturbation equations. This 
situation is very similar to the usual vacuum energy (cos- 
mological constant) that gravitates (and therefore con- 
tributes to the Hubble equation) but does not contribute 
to the perturbation equation. Another example of this 
kind is the Casimir energy, that represents only a global 
shift in the vacuum energy, which can be negative and 
therefore resolve the singularity without effecting per- 
turbations [13, HI]. Yet another interesting scenario is 
to consider higher derivative gravity correction terms to 
the Einstein action which we expect will become dom- 
inant near the bounce, as long as the curvature at the 
bounce point becomes comparable to Planck-scale cur- 
vature. Recently, a ghost-free and asymptotically free 
higher derivative gravity model leading to a cosmological 
bounce was proposed in Q , and in a follow-up paper ^] 
it was suggested that a long bouncing phase may lead to 
the correct thermal string gas initial conditions for the 
new structure formation scenario [3l|. [s^ l (see also [s^ 
for a recent review). The action of [5| has the property 
that, although the cosmological background is changed, 
the correction terms in the equation of motion for long- 
wavelength (compared to the scale of new physics) cos- 
mological perturbations are suppressed, and we can thus 
evolve these perturbations using the usual equations for 
cosmological fluctuations ^ . 

Thus in this paper, we will study the evolution of 
cosmological fluctuations using the perturbed Einstein 
equations through a specially chosen cosmological bounce 
background which leads to exactly soluble equations. In 
our example, we find that the growing mode in the con- 
tracting phase goes over into the growing mode in the 
expanding phase, in contrast to what happens in the 
model of [18|. Our result gives rise to the hope that, 
once the cosmological singularity in the Ekpyrotic sce- 



Although such a bouncing cosmology can provide us with a non- 
singular description of the universe and possibly a new mech- 
anism for the origin of cosmological perturbations, the models 
have not yet addressed the question of why after the bounce, in 
the expanding branch, we are mostly left with Standard Model 
particles. One has to perhaps find a mechanism similar to "re- 
heating" where after the bounce most of the energy is converted 
to Standard Model degrees of freedom. Another possibility could 
be that the bounce creates the initial conditions to drive a phase 
of late time inflation [53 . l35l | and/or create Standard Model 
baryons and cold dark matter by exciting the minimal super- 
symmetric degrees of freedom j3a | (for a review see [37j'). 
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nario is resolved by a nonsingular bounce, it will not be 
necessary to invoke a scale- invariant spectrum of initial 
entropy fluctuations [H, [M SO, SH SI S, El in addi- 
tion to the existing scale-invariant spectrum of adiabatic 
modes in order to obtain a scale-invariant spectrum in 
the expanding phase. 



Thus, for the purpose of illustration and technical sim- 
plicity we choose to work with a specific ansatz^ to cap- 
ture the essential effects on the perturbations as they 
evolve across the bounce: 



n 



qr 



(7) 



II. THE MODEL 

The starting point of our analysis is the generalization 
of Einstein's equations to 



(2) 



where Qf^u arises from some uknown new physics respon- 
sible for resolving the Big Bang singularity. A priori. Qui/ 
may arise either in the gravity sector (see for instance [5|) 
or in the matter sector (26l l27l [H, [2^ but we purposely 
do not specify it. As usual the above equation can be bro- 
ken down into the background (which only depends on 
time) and perturbed (depends on both space and time) 
quantities such as 



(3) 



and similarly for and QJ^. For homogeneous isotropic 
cosmology, the background quantities generalize the 
Hubble equation to 



3i?^ =< Gl >=< To" > + < g° >=< p> + < 



(4) 

To keep things as general as possible let us not specify 
the new physics, i.e. < Qq > that causes the universe 
to bounce, but simply make an ansatz^ for the time evo- 
lution of the Hubble rate. Near the bounce, the most 
generic behavior of the scale factor is given by 



(5) 



a(T) = 1 + - 



where tq corresponds to the bounce time scale. ([5]) leads 
to a Hubble evolution of the form 



n 



(6) 



We note that this has the right asymptotic property as 
T — > ±cxD, while during the bounce 7i goes linearly with 
T as expected. 

As matter we consider a fluid with an equation of state 



P 



ujp, 



(8) 



where p and p stand for the energy density and the pres- 
sure of the fluid, respectively. Inserting this equation of 
state into the general relativistic Friedmann equations 
yields the following time evolution of the scale factor: 



air) = aolrj^ with q . 



l + 3aj 



(9) 



where ao is a constant. 



III. PERTURBATIONS ACROSS THE BOUNCE 

Starting from ^ and using perturbative expansions 
such as ([3]) we obtain the generalized perturbation equa- 
tion: 



sg; = + 6q; 



(10) 



We now make the crucial assumption that SQ'j^ vanishes 
(or remains negligible during the bounce). As mentioned 
in the Introduction, when the resolution of the Big Bang 
sin gula rity involves new non-local physics such as @, |26l 
[27l. l28l. [soj . then at least for fluctuations whose length 
scales during the bounce are much larger than tq, this 
can be justifled. 

Returning to pil)) . one can now check that just by us- 
ing the usual expressions of 6G'j[ and (5T^ in terms of $ 
and Sp for an ideal gas, one can derive the usual Gen- 
eral Relativistic perturbation equation for the Bardeen 



where Ti. is the Hubble rate in conformal time. Now, 
we also know that at late times H has to have the right 
asymptotic property, 7i — > with q < 1. Of course, 
in practice the precise transition from the accelerating 
bouncing region ([6]) to the deccelerating late time regime 
will depend on the new physics that is introduced. How- 
ever, at least in the case where there is really only one 
new fundamental scale we expect only the details to de- 
pend on the precise functional form of the transition. 



If one wants, one can deduce < QjJ > from tlie ansatz. 



^ If we believe that the singularity at r = will be resolved, then it 
is clear that the pole at r = will have to shift from the real axis 
to the complex plane. Since we want Ti. to be finite and real in 
the entire real axis, it is easy to see that in fact the poles must lie 
along the imaginary axis and come in conjugate pairs. Therefore, 
the simplest ansatz that one can make for H, is to consider only 
a pair of simple poles as in JTJl. However, the general algorithm 
that will be advocated here can be carried forward for any non- 
singular bounce ansatz that one may get for a given new physics. 
Also, we believe that this ansatz for the Hubble parameter may 
be realized as non perturbative solutions of non metric theories of 
gravity as well as from four dimensional, higher order corrected 
Heterotic M theory actions. These constructions will be pursued 
in future work [29| . 
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potential $ 

$fe+3(l+cj)7^$fc+cjfc2$fc + [27^' + (l + 3cj)7^2]$A- = 0. 

(11) 

Crucially we note that this derivation does not "re- 
use" the background Hubble equation and therefore goes 
through even though the background equation is mod- 
ified. As we have expalained in the appendix, this is 
not the case for the standard evolution equation for the 
Mukhanov variable (also see section ??) and in fact this is 
the reason why one cannot trust the General Relativistic 
equation for the Mukhanov variable to track perturba- 
tions around the bounce. 

In passing we note that for GR, in case of an ideal fluid 
the last term drops out. However, when one modifies 
gravity, this is no longer true, and this is what makes the 
analysis a lot more interesting. 

Substituting our ansatz ([7]) into pip , we obtain the 
following equation for the metric fluctuations 



3{l+Lu)qT 



$fc = 0. (12) 



As one can see, the coefficients in the above differen- 
tial equation are all non-singular and completely well 
defined. Therefore one can in principle solve for in 
the interval (— cx), cx)) and understand how an initial per- 
turbation evolves from the contracting to the expanding 
phase, without having to implement any matching condi- 
tions anywhere. Although this is in general only possible 
to achieve numerically, as we will argue now, one can 
make significant progress analytically. 

Let us first identify some of the regimes where the dif- 
ferential equation may simplify. Firstly, one has the usual 
division between the super- and sub-Hubble regimes sep- 
arated at points where 



k = \H\ ^ k\T\^ -q\T\+T^k = 0. 



(13) 



As is typical in bouncing cosmologies, there are four so- 
lutions to the above equation: 



± 



4fcro 



2fc 



±T±. 



(14) 



Two occur during contraction (— t±), and two during ex- 
pansion (tj-). The initial and final crossings are the usual 
ones for a contracting and an expanding universe, re- 
spectively, whereas the two middle ones occur near the 
bounce point and are a special feature of nonsingular 
bouncing cosmologies. The first crossing occurs at — t+, 
when GR is a valid description of the background and 
modes exit the Hubble radius (note the universe starts 
from a cold 7i ^ phase, so that all the modes are sub- 
Hubble to begin with). Around r ~ tq, the new physics 
kicks in and 7i starts to decrease again. Thus, all the 
modes which are super-Hubble start to come back inside 
the Hubble radius, and for a given mode this happens 
at — T_. By the time the bounce occurs all the modes 



are again sub-Hubble, but after the bounce as the Hub- 
ble radius reaches a finite value, modes again exit the 
Hubble radius. This happens at t_. Finally, as the ex- 
pansion continues according to the equations of standard 
cosmology, modes eventually enter our horizon at tj^. 

A second way to divide the evolution is to consider the 
range |t| > tq and |r| < tq separately, where the former 
corresponds to considering the intervals when GR is a 
good description of the background, while for the latter 
one has to include the effects of the new physics. So, how 
should one proceed? 

First let us look at as |r| ^ tq, where we recover 
the general relativistic limit 



" , 6(1 -Hw) ^, , 



$ , + $fc = , 



(15) 



r(l + 3tj) ^'^ 

which has the following analytical solution 

$fe = T-'^[A:-''i?±(fc)J„(V^fcr) + r5±(fc)J_,(v^fcT)], 

(16) 

where 



1 /5 + 3t^ 



2 V 1 + 3w 



(17) 



and -I- or — labels the coefficients corresponding to the 
expanding (tq, oo) and contracting (— oo, — tq) phases re- 
spectively where the solution is a good approximation. In 
particular, it incorporates the sub-Hubble phases at early 
and late times, i.e. between (— oo,— t+) and (r-|_,oo). 
Also, we have included the fc-dependent factors in front 
of the coefficients D± and S± to ensure that in the super- 
Hubble phase the k dependence is completely contained 
in Diki) and S[ki). By using the asymptotic form of Bessel 
functions: 



lim Jvix) 



(18) 



and taking the appropriate t — > limit of p6p . this can 
easily be seen: 

lim $fe = D±{k)u''^^ + 5±(fc)cj-"/2(±r)-2- . (19) 



Now of course, if we know the solution to in the 
entire range (—00, 00), then we will know how {-D+, 5*+} 
is related to {-D-, "S*-}, but as we will now explain, one 
requires less. Finding this "transfer matrix" without any 
ambiguities is in fact the main endeavor of this paper. 

In order to achieve this let us next look at the super- 
Hubble phases. As usual, in these phases, the fc-terms 
can be ignored, asTi. > k and the second term dominates. 
Finally, we are left with the sub-Hubble phase around the 
bounce and now comes a crucial observation. Unlike in 
the usual GR scenario, in this phase which occurs once 
the new physics has kicked in, i.e. \t\ < tq, the /c-term 
can again be ignored in favor of the fourth term. This is 
because for the modes that we observe today, typically^ 



If there is no prolonged inflationary phase. 
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k <C Tq"^, the latter being the scale of new physics which 
is expected to be close to the string or the Planck scale. 
As a result, we find that for the entire range (— T+,r+) 
we can ignore the fc-term and need only solve 



3(1 + uj)qT^, 



(r2 



^fe = 0. 



(20) 



A few comments are now in order. Firstly, observe that 
we have completely bypassed having to do any match- 
ing at ±T_ or at r = 0, as is often done in literature. 
Secondly, the formalism described here divides the en- 
tire evolution into three overlapping regions: (— oo, — tq), 
(— r_|_,r_(-) and (to,oo). As a result, the matching of the 
solutions to (fT^l) and ([^0)1 can be done unambiguously 
by looking at the asymptotic properties of the solutions. 
This is to be contrasted with matching conditions which 
are imposed at specific points in time (like between sub- 
and super-Hubble phases at the point of Hubble cross- 
ing). A final remark, as is clear the sub-Hubble phase 
around the bounce is very different from the usual sub- 
Hubble phases (for instance in inflationary cosmology) 
because the fc-term is unimportant, and thus any analy- 
sis based on the usual framework of solving the equation 
keeping only the fc-term is bound to give incorrect results. 

Above, we already obtained exact solutions for the 
1st and Ilird region, i.e. in the intervals (— cx), — tq) and 
(ro, oo). In order to avoid any ambiguity involving mode 
matching we need to find exact solutions to ([^01) . To do 
this it is convenient to define 

y.^Mr'+T^)-. (21) 

The differential equation then becomes 

(t2 + T^)^y" + (r^ + T^)T{An + 3(1 + Lu)q)y' 
+ [2n{T^ + T^) + {4n{n - 1) (22) 
+6nq{l + uj))t^ + 2qTl]y = . 

If we now choose n to satisfy 

An{n - 1) + &nq{l + uj) = 2q , (23) 

which leads to 

-1± V2 + 3w 



n± = 



(24) 



and in particular choose to work with n — n_, then the 
ordinary differential equation greatly simplifies and one 
has 



IV. MODE SWITCH: AN EXACT SPECIAL 
CASE 

In order to have a clear understanding of the physics 
we will focus on a special case where the solution can be 
written in terms of familiar functions. This is the case 
when the coefficient of the second term, /3, equals one, 
which happens when uj « 5.314. 

We note in passing that such a large u can actually be 
physically interesting as it is known that as w ^ oo we 
produce a scale-invariant spectrum in one of the modes 
(the Ekpyrotic scenario). In any case, the reason we want 
to consider the above uj is because it lends technical sim- 
plicity to extract the physics which will be useful when 
we study the more general case. 

For the above special case the differential equation sim- 
plifies to 



(t^ + T^)y" + Ty' - a^y = 



with 



2ni > 0. 



(27) 



(28) 



Its solution reads 



y = Bi exp ( a sinh ^ — ) + B2 exp ( —a sinh ^ — 
V Taj \ To 

(29) 

In order to understand how the perturbations propagate 
through the bounce we have to look at the behavior of y 
as T — > ±00. Now, 



lim sinh 

r — ^±00 



±ln 



±2r 



,T0/ V Tq 

As a result, the asymptotic values of y reads 

'±2t\±" /±2t\^" 
+ B2 



\im y = Bi 

T— >±00 \ To 



V ^0 



(30) 



(31) 



Something remarkable has happened: in going from 
contraction to expansion, the coefficients corresponding 
to the two modes r" and t~" have completely switched. 
The dominant mode in the contracting phase goes over 
into the dominant mode in the expanding phase, un- 
like what happens in a singular bounce making use of 
the Hwang- Vishniac and Deruelle-Mukhanov matching 
conditions. To see this more precisely, by matching the 
bouncing solution (|29p to the late time solutions given 
by (fT6| in the overlapping regions one finds the exact 
relation: 



where 



{T^ + T^)y" + T(3y'-2n+y = 0, (25) 



2(1 + 3uj- V2 + 3w) 
l + 3w 



(26) 



and S+ = D_ 



(32) 



For a general equation of state, although we do not 
expect such a complete switch, the above result certainly 
suggests that there would be mode mixing. 
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V. MODE MIXING AS a; > 1 

It is known that the Ekpyrotic scalar which mimics 
an ideal fluid with a large equation of state parameter 
produces a (nearly) scale-invariant spectrum of pertur- 
bations during the phase of contraction, but in the grow- 
ing mode. It was argued that this mode matches to the 
decaying mode exclusively during the expansion phase 
and therefore cannot explain the near scale-invariance 
observed in the CMB today. The previous exact analysis 
already suggests that this may not be true if the Ekpy- 
rotic bounce is smoothed out. 

Quite remarkably, in the a; — > oo limit the differential 
equation simplifies to 



{T^+T^y" + 2Ty' = 0, 



(33) 



and again becomes amenable to an exact treatment. The 
equation has the following solution 



y = Bi + Satan- 1(t/to). 
Now, in the case w ^ 1 we obtain from ^9, 

Vk 



$fc = D±{k)uj^^'^ + S±{k)uj-^^^{±Ty\ (35) 



where we have used that n —^ as to oo. We can now 
relate the -I- coefficients with — via the bounce solution 
By considering the r ^ ±oo limit of one easily 

finds 



u;^/^D± = Bi± 



This leads us to 



and ±uj'^/^S± ^ ~B2Tq . (36) 



D+(k) = D(k) -j^S-{k) and S+{k) = -5_(fc) . 

(37) 

We notice that the constant mode -D+ gets contributions 
from both the modes in the contracting phase. In par- 
ticular, the contribution from S- gives a scale invariant 
spectrum ^! 



VI. DISCUSSION 

In this paper, we have followed the evolution of cos- 
mological fluctuations across the bounce using the rela- 
tivistic potential <&, the metric fluctuation in longitudi- 
nal gauge. From experience built up in investigations of 
inflationary universe models (see e.g. [lE[3])j Pre-Big- 
Bang cosmologies and the Ekpyrotic scenario, working 



^ To see this note that according to conventional quantization of 
the Mukhanov variable, as r — > — oo, v ~ {Vb/v^)[cos(v^fcr) + 
sin(^A:T)]. Using the relation ~ t3(i+")/(i+3'^) [$'+2!/($/t)] 
valid for ideal fluids in the General Relativistic limit (— oo, — t+), 
one then finds D— ~ k~^/'^ and S— ~ k~^/'^, the latter giving 
rise to a scale-invariant contribution to the power spectrum. 



in terms of can be dangerous because of very sensi- 
tive dependence of the evolution on matching conditions 
at times when the equation of state of the cosmological 
background changes, e.g. during the period of reheating 
in the context of inflationary cosmology. In our example 
we are safe from this danger because we have an exact 
solution and hence do not need to invoke any matching 
conditions. 

On the other hand, from the point of view of the quan- 
tum theory of cosmological perturbations (see [47,, ,4^ 
for pioneering works), a variable different from $ carries 
more physical meaning, namely the variable v in terms of 
which the action for cosmological fluctuations has canon- 
ical kinetic term. In the case of hydrodynamical matter, 
the variable which determines the curvature perturba- 
tion in comoving gauge, is related to $ via 



c,v = u 



(34) where Cg is the speed of sound. 



and 



a \ p 



-1/2 



VP + P 

The variable v obeys the following fleld equation: 



v" -c^V^v w = 0, 

z 



where 



-ncs 



and p EE n^-n' . 



(38) 



(39) 



(40) 



(41) 



(42) 



In inflationary cosmology, the variable v remains con- 
stant during the transition of the equation of state which 
takes place at the time of reheating, whereas $ jumps by 
a large factor. It is thus the variable v which is a more 
robust one to follow. In fact, in the absence of entropy 
fluctuations, one can show that on super-Hubble scales, 
the variable 



7^ = - 

z 



is conserved |45l Wo, k 



(1 + w)7t 



0. 



(43) 



(44) 



In the contracting phase of an Einstein universe it can be 
shown [l2| that the dominant mode of $ does not couple 
to the variable v, and hence the curvature fluctuation is 
only sensitive to the decaying mode of a mode which 
does not lead to a scale-invariant spectrum. 

In a bouncing cosmology, it becomes problematic to 
use the equation of motion (^1]) for v. For instance, at 



7 



the bounce point the variable z blows up and hence the 
equation is singular. Essentially the new physics which 
solves the Big Bang singularity necessarily has to modify 
(|4T|) . As we discuss in details in the appendix, in general 
it is not even possible to come up with an equation for v 
(see [i^l for an attempt in this direction) while at least 
in our scenario the perturbation equation for $ is still 
correct. Thus it is not justified to use the large scale 
analysis and conclude that 7?. is constant. These are our 
reasons for focusing on the evolution equation for $ and 
not for V. 



VII. CONCLUSIONS 

In this paper we have presented a particular nonsingu- 
lar bouncing cosmological background, in which the usual 
general relativistic equations for cosmological perturba- 
tions can be solved exactly. We find that the growing 
mode in the contracting phase goes over into the domi- 
nant mode in the expanding phase. This is unlike what 
happens in four dimensional effective field theories of Pre- 
Big-Bang or Ekpyrotic type, in which the fluctuations are 
matched at a distinguished but singular hypersurface us- 
ing the analog of the Israel matching conditions. Our 
result supports the conclusions of [ll,|l7| that the trans- 
fer of fluctuations is very sensitive to the details of the 
bounce. 

Our analysis assumes that terms in the Lagrangian dif- 
ferent from the Einstein-Hilbert terms generate the non- 
singular bounce. At the same time, it is crucial that these 
new terms not effect the equations for the IR fluctuation 
modes. An example where both conditions are satisfied 
is given in [H, |33| ■ 

Applied to the Ekpyrotic scenario, our result implies 
that one may not be required to invoke entropy fluctua- 
tions with a scale- invariant spectrum [H, [s^ [40, IHI, HI] 
in order to obtain a scale-invariant spectrum of curvature 
fluctuations in the expanding phase. 
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VIII. APPENDIX 

In this appendix we explain why in our bouncing uni- 
verse scenario we used the Bardeen potential $ to track 
the perturbations across the bounce rather than the more 
conventional Mukhanov variable v. The discussion will 



also clarify how we can evade the usual conservation ar- 
guements involving the curvature perturbation (|44p to 
get mode-mixing. Let us first define the intermediate u 
variable: 



u = exp 



1 + cDn dT 



(45) 



For an ideal gas using the continuity equation 

p' + 3H{p + p) = (46) 

u simplifies to 



If one now uses the relation 



27i:' + (I + 3cj)7i;2 = 



(47) 



(48) 



which is valid only for a General Relativistic background 
one can derive a relatively simple differential equation for 
u from pT|) : 



- ujVu - 



= where 9 = 



1 



a\/l + LU 



(49) 



Since the new physics that resolves the singularity al- 
ters the background equations, (|48l) can no longer be 
valid (in fact one can check that (I48|) always lead to a 
singular universe), and therefore one cannot trust (j49p . 
Nevertheless, one can actually derive a similar equation 
for u which does not use and therefore will is valid 
irrespective of the new physics. The equation reads: 



wv^u- {an' + pn^)u ^ 



(50) 



where 



a.^and/3.^^!±^ (51) 



One can easily verify that provided is valid, ([50]) 
reduces to (02]). 

In the usual General Relativistic analysis one typically 
defines the Mukhanov variable via ((38| which for an ideal 
gas simplifies to 



V ^ u' + Hu 



(52) 



and turns out to be the right canonical variable for the 
effective quantum action [ol Il0|| (again in the General Rel- 
ativistic limit) and therefore let's one impose the "sub- 
Hubble quantum fluctuations" as initial conditions. In 
this context we note that using the relations between 
<&, u and V (|40l38p . one can also read off the appropriate 
initial conditions for $ or m, from that of v, and then 
choose to analyse the propagation of the fluctuation in 
any of the variables that one chooses to (see footnote 
9 for instance). We should also point out that all the 
definitions (|40l38p are non-singular, and one expects the 
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solutions for all these variables to be regular when one 
uses the correct evolution equations for them. 

Actually, this is where the Mukhanov variable as de- 
fined in ([38)1 looses it's usefulness. Unhke (j) and u, one 
cannot even find an evolution equation for v like (|4ip 
in the general case when one has introduced some new 
physics to resolve the singularity. Instead one finds the 
following equation 

v" ~~ toV^v - [(a + 2)n' + pn^] {u' + fu) = (53) 

where we have defined a function which depends on the 
background: 

^ ^ (3H^ + (q + l)n" + (2(3 + a)HH' ^^^^ 

It may be an interestin whether it is possible 

to generahze the definition of v such that it reduces to II41I I in the 
GR Hmit, but differs from it near the bounce in such a way that 
one can obtain an evolution equation for v like it's GR cousin. 



Clearly, in order to have an equation only in terms of v, 
the terms involving u must combine to give us v which 
can happen iff one can satisfy 



3{1 + uj)H" + [1 + 3u} fnn' ^ (55) 



Indeed for GR this condition is satisfied. However, it is 
clear that when one introduces new physics, in general 
there is no simple evolution equation for v and (j41|) is 
certainly not valid^". It also explains why one can't use 
the conservation law (144]) which is based on (|4T|) . This is 
why a mode-mixing is actually possible in our scenario. 

For such an attempt in the context of quantum cosmological 
models see [53 |. 
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